NACA TN 2447 



NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 


TECHNICAL NOTE 2447 


} ^ d - ***- tos * 

STATISTICAL FLUCTUATION OF INTENSITY IN DEBYE -SC HE ERE R 
LINES DUE TO RANDOM ORIENTATION OF CRYSTAL GRAINS 

By Hans Ekstein 
Armour Research Foundation 


v v 

NACA 


Washington 
August 1951 


mm 

fU&jUf.,.? ] i 5 

ill 5 Sf ? 4- • £ 


l ummi 


TECH LIBRARY KAFB, NBA 



I Cvn UDnHni rwre, nhi 



DQbSSM? 

NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE 2447 


STATISTICAL FLUCTUATION OF INTENSITY IN DEBYE- SCHERRER 


LINES DUE TO RANDOM ORIENTATION OF CRYSTAL GRAINS 


By Hans Ekstein 


SUMMARY 


In the Debye-Scherrer diagram of a stationary polycrystalline 
sample, the intensity distribution in the line is erratic because of 
statistical irregularities of the crystal grain orientation. The fluc- 
tuations of intensity are discussed for the case that the only causes 
of line breadth are the natural spectral width of the primary radiation 
and the finite size of the grains (i.e. , pure Fraunhofer diffraction). 
The mean deviation of the centroid of the intensity curve from the Bragg 
angle for the strongest wave length can be considered as the error in 
measuring the line. This deviation is theoretically estimated. For a 

-4 

typical case, namely, grains of 2 x 10 centimeter length and 1-square- 
millimeter irradiated area, the error is found to be about 10 times 
smaller than that reported in experimental literature. 


INTRODUCTION 


A stationary polycrystalline aggregate of crystals gives rise to 
an intensity distribution in the Debye-Scherrer line which depends on the 
number of grains being in reflecting position. In general, a line of 
erratic shape is obtained with no well-defined maximum. This difficulty 
is usually avoided either by rotating or oscillating the sample about 
an axis normal to the incident beam (thus averaging out the statistical 
irregularities) or by choosing the effective number of grains so large 
that the relative fluctuations become negligible. 

In some cases, however, neither of these devices is practicable. In 
particular, for stress measurements, sample rotation (or, what amounts 
to the same, oscillation of the tube) is often cumbersome, and the grain 
size of the sample is not subject to choice. But in other cases too, for 
instance, with focusing arrangements, the irregularities of crystal 
orientation cause difficulties, particularly when high resolution is 
obtained by the use of high diffraction angles or large distances between 
sample and film. 
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In the following, the error in the measurement of the line center, 
due to these irregularities, is estimated. For this purpose, a definite 
geometric condition must he chosen which determines the intensity dif- 
fracted by a single crystal grain in a general orientation. In most 
cases, the width and divergence of the incident beam will be the main 
determining factors for this intensity. However, for very high resolu- 
tion, these geometric factors become small in comparison with the width 
due to the spectral impurity of the primary radiation (reference l) and 
it is in this case precisely that the question of accuracy becomes most 
important. Consequently, a primary beam of negligible width and diver- 
gence is assumed so that the recorded intensity depends only on the angle 
of diffraction; that is, pure Fraunhofer diffraction is considered. The 
crystal grains are assumed small enough so that the kinematic theory can 
be used. The crystal shape is' assumed spherical with a uniform diameter. 

The investigation of the method reported herein was made at the 
Armour Research Foundation under the sponsorship and with the financial 
assistance of the Rational Advisory Committee for Aeronautics. This 
report is part of a cooperative project with Dr. Stanley Siegel. 


DIFFRACTION BY A SINGLE CRYSTAL GRAIN 


In figure 1, 0 is the origin of the reciprocal lattice, and B 

is the position of the reciprocal lattice point belonging to a Bragg 
lattice distance d. For a given crystal lattice, B can have any 
position on a sphere of radius l/d and center 0. The Laue spheres 
belonging to different wave numbers k^, 1 q, and k£ have the centers 

Ff, Fq, and F2> respectively. If Rq is the mean wave number of the 

primary radiation (corresponding to the maximum spectral intensity), 
then the strongest reflection will take place when B is situated on 
the intersection of the Laue sphere with center Fq and the "reflection 
sphere" with center 0 and radius l/d. The diffraction maximum will 

occur at the angle d Q so that [_ HOC = ~ - t 3 q. • The intensity will 

decrease rapidly as d deviates from this value. Similarly, positions 
of B which are far from the intersection (point C of fig. l) will 
contribute little intensity. Therefore, all angular distances from the 
point C can be considered small. 

Returning to a general wave number k, consider the intensity of 
wave number k diffracted in a direction given by the point A and 
caused by a crystal position indicated by the point B." Let B have the 
polar coordinates l/d, it/2 - d + 7, and 0 with respect to the origin 
0, and let the point A have the polar coordinates 2k sin d, jt/2 - d, 
and Tie. With the assumption of spherical shape, the intensity will be 
a function of the distance R = AB. 
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The following equation is obtained 


cos (-3-7) 

( 1 ) 

2 

If \jr is the interference function, the intensity of wave number k 
diffracted in the direction 8 is 


R 2 = 


+ 4k^ sin 2 8 


- £ 


sin 8 cos: 7 + (l - cos \|r) sin 28 


J k - 5 2 (R) (2) 

Let the spectral intensity distribution be given by the function f (k) . 
Then the total intensity of all wave lengths is proportional to 

J* f(k)t 2 (R) dk (3) 

where R is considered as a function of k according to equation (l). 
In most measurements of this type, the film is rotated about the primary 
beam so that intensities belonging to a definite -8 are averaged with 
respect to the azimuthal angle. This fact may be taken into account by 
integrating the expression (3) over all azimuthal angles of the 
reciprocal lattice point. The intensity becomes: 


J 



d t dk f(k)$ 2 (R) 


(4) 


The exact form of the interference function ijr was calculated by 
Patterson (reference 2) . It shall, as usual, be replaced for math- 
ematical convenience by a Gaussian function Ce" 0 * 2 of identical ratio 
of maximum to the integral of the squared function. This requires 
(reference 2) that 

= p2 *0f) = 8 ' l6p2 (5) 

where p is the radius of the crystal grain. For mathematical con- 
venience, the resonance curve f(k) which .represents the spectral s hap e 
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of "the characteristic radiation will also "be replaced "by a Gaussian 
curve so that the two curves have the same half-value width. If w 
the measured half-value width in the scale of the wave lengths 
function f(k) will he 


f(k) 


( k - k o f + (V^ 2 ) 2 


and this function has the same half-value width as 


e -(3(k-k 0 ) 2 


if 


k lp ge 2 = 2.7 6 
w 2 ^^" 


Equation (4) now takes the form: 


J 




.. .. -a R 2 -p(k-ko) 

ch|r dke e ' ' 


2 




-a(l-cosijf) sin2-acos(-fl-7) 
e 


dijr 



I 0 (agk)e agk e 


•p(k-k 0 ) 2 -ako S 


dk 


where 


g = 


2 

d 


cos - 7) sin 2£ 



+ 4k 2 sin 2 Q 
dr 


4k 

d 


"Tti. * ft 

-T~ sm d 


cos 7 


is 

the 


( 6 ) 


(T) 


(8) 


(9) 


(10) 
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Since only the close surrounding of k = k Q will contribute substantially 
to the integral (8), the limit of integration may be changed to 

Furthermore , if sin 2d is not too small, the function I 0 (agk)e" a S k 

is slowly variable and can be replaced by the constant I Q ^agk-^e -0 ^! 
where 

2ct 

PkQ + -jj- sin d cos 7 

k l p (11) 

3 + a sin d 

is the value of k where the exponential factor in the integrand attains 
its maximum value. 

There can be written 





dk 


( 12 ) 


or by equation (10) : 


J oc I 




* e- D 


P + 4x sin 2 d 


( 13 ) 


where 


D = pk Q ^ + 3u 


a ^2gkQ + lj£ sin d cos 7^ 
4^P + tac sin 2 dj 


( 14 ) 


Putting 


sin d = 


2k 0 d 


( 15 ) 
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sin d = sin t9q + A sin d (i 6 ) 

and 

cos 7 = 1 - ^ 5 - (17) 

expression (l4) is expanded up to second-order quantities (A sin $) 2 
and 7 2 to obtain: 


D = tack 2 / 7 2 sin 2 d Q + \ (l 8 ) 

l tat. sin 2 d 0 + 3 / 


SHAPE OF DIFFRACTION LINE 


The spatial orientation of a crystal grain is characterized by the 
position of a crystallographic direction on the unit sphere. If the 
direction (hkl) corresponding to the reciprocal lattice point B is 
chosen as characteristic for the orientation of a crystal grain, then 
the set of representative points on the unit sphere will determine the 
diffracted intensity. The unit sphere is subdivided into small cells of 

area A 2 (fig. 2). 

The intensity diffracted by directions (hkl) which form angles from 
cp to cp + A with the primary beam will be proportional to the number 
of representative points situated on a circular ring of width A 
(fig. 2). If this number is n^', the total intensity diffracted in a 
direction $ by planes (hkl) is 

21 


where J is the intensity given by equation (13) • According to figures 1 
and 2 , the connection between 7 and cp is 

cp = 7 + — - d 


<4 


(19) 
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Since there are always crystallographic directions which, for diffraction, 

are equivalent to a given direction (hkl) , ^ all contributions to the 
intensity have not yet been counted. The contributions of all those 
crystals whose equivalent directions form an angle cp with the primary 
beam should be added. Since only random orientations have to be dealt 
with, this is done easily by multiplying the n umb er n^* by s, the 
number of equivalent planes. If 


= V ‘ S ( 20 ) 

the total intensity becomes 

J t = (21) 

9 

The exact shape of the intensity curve is unimportant, but determining 
that angle -d which experimentally appears to be the center of the line 
is necessary. This requires the definition of the experimental procedure 
by which the center of an irregular diffraction line is measured. In 
fact, most measurements of this type are visual. For the present 
purpose, a more quantitative method, at least in principle, is obvi- 
ously needed. Actually, it may well turn out that the more refined 
methods of objective measurements are unnecessary because the inherent 
error is large. However, the aim at present is to calculate the best 
possible accuracy with a given crystal aggregate. 

Probably the best experimental procedure would be to find, by the 
method of least squares, that resonance or Gaussian curve which fits 
best a measured intensity curve, but this method is so cumbersome that 
the line center has been defined herein as its center of gravity or 
centroid | which is given by 



( 22 ) 


The integration should be extended over "the line" an ill-defined 
quantity since the intensity between two lines is the sum of the inten- 
sities belonging to the two, or even more, nearest lines. This difficulty 
is typical for the precise characterization of a line. The definition (22) 
is meaningful only if the intensity between two lines decreases to 


According to Friedel's law, there is at least one equivalent 
direction (-h, -k, -1). 
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a negligible value, in which case the integral must he extended into 
this region of substantially zero intensity, and the exact position of 
the limiting point does not matter. It must be assumed that such a 
region exists. 


CRYSTAL STATISTICS 


The N crystal grains of identical shape are oriented at random. 
The probability that there are ru grains whose directions (hkl) form 
an angle qp to cp + A with the primary beam is that of n^ balls out 
of N falling into a box of area 2 it A . sin cp when they are thrown at 
random into a total area hit. It is 



— sin cp 
2 


N-n<p 


N! 

n^KN - n^TT 


( 23 ) 


As explained above, the existence of s equivalent lattice planes 
modifies the situation by substituting the effective number sW of 
grains for N. The mean number njp is 



K sin cp 
2 


( 24 ) 


and, from equation (23), one concludes in the usual way. 


M 2 - - “(p) 2 


— sA ' N sin cp 
n cp " " q 2 


( 25 ) 


The mean square deviation of the position of the centroid from its' 
mean value is now considered. If is taken as origin of the angular 
coordinate, the mean value of £ is zero, and 


t 



cl(Afl) 

O 



J t (fl) d(Afl) 


(26) 
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The coefficient of g is the "integrated intensity" of the line. It 
is also subject to statistical fluctuations hut clearly these fluctua- 


tions are of minor importance for the value of g 2 , and A d(Ad) may 
be replaced by its mean value. The following expression is obtained for 

O 

g by equation (21): 



d(Atf) 2 | 2 


(A-8)J t U) d(Ai3) 


w 


J(<Pitf)(Ad) d(Ad) 


( 2 ?) 


where the index i corresponding to the angle cp^ has been substituted, 
for the general index cp. 


If 


J(q>i,d) (Ai 9) d(Ad) = f(cpi) 


(28) 


and 


then 



(29) 


cf = n-^ffcp^ffcpT,) = ^nin^cpj^cpt) ( 30 ) 

i,k ' i,k 


In the double sum (30), the terms with i = k and those with i ^ k 
are separated so that 

Ci 2 = YL ^ i n k f(q3 i )f(cpk) +Xn i 2 f 2 (cp 1 ) 
i A i 


and, since by hypothesis the numbers n^ and n k for i / k are 
entirely independent, 
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C|‘ 


2 = 51 14 n k f(<pi)f(<Pk) +2Z n ± 2 f 2 (cp i ) 




(31) 


It is remembered that 


n i 2 = ( n i ) 2 + (^i ) 5 


so that, by equations (24) and (25 ) , 

<2. 


C£ 2 = y * sin cp ± sin cp k (A ♦ N) 2 f (<Pi)f (<Pk) 

ij^k ' * 


o. ./ sA * N sin cp,-\ sA * N sin q>i 

f 2 ( 9l) i + — ij 


(32) 


It is now possible to pass to the limit A — >0: 


Cl 2 - 


• (?)'(/ 


?/• 


sin cpf (<p) dcp ) + / sin cpf 2 (9) dc P 


(33) 


since that part of the second sum in equation (32) which carries the 
factor A 2 vanishes in the limit. 

Further simplifications are necessary in order to derive a simple 
result from equation (33). In equation (13), the factor preceding the 

Gaussian term e _D will be considered constant since it is slowly 
variable in comparison with e -D . Equation (13) becomes: 


-4aV 


J « Ke 


p - 2 P(Asin^) 2 

7 3in *° 4asin 2 V3 


(34) 
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If 


Ap = cp - cp Q = cp 



by equation (19) , 7 becomes 


7 = A3 + Ap 

Furthermore, if 3 is not too large, 

« 

A sin 3 ~ cos 3 qA3 

Temporarily, put 

iktk Q 2 = A 

. p 

sir 3 q = a 
P cos 2 3 n 

_ — = C 

4x sin 2 3 q + P 

Ap = x - 
A3 = y 

so that equation (13) can be -written 








-A 

J «Ke 


ja(x+y) 2 +cy^j 


(35) 


(36) 


(37) 


(38) 


(39) 



12 


WACA TH 2bk r J 


By equation (28) 


f(x) = K r e.‘ A E (x+y)2+Cy ?ly dy = - -JEL. lZZjZZ 
J.oo c + a 1/A(c + a) 


Aacx d 

c+a (to) 


Since f is an odd function of Ap, the first integral in equation (33) 
vanishes: 


2 _ Ns K 2 ; 


. 2A.ac(Ap)‘ 


cr = 


A(c + a) 


2 P 00 o • 

EJl — / sin qp * (Ap) e c+a d(Ap) 

+ a) ^ 


(to) 


Because only crystals with angles <p rj cPq, that is, small Acp, will 
contribute to the intensity, sin qp % cos $q may he considered as a 
constant in the integral of equation ( to) : 


C £ 2 = -j[Egf.qL_ cos l/fl /c + aW 2 
2A(c + a) 3 ° 4 \2acA / 


KtTs cos 

ft 

OJ 

CO 

2ac(c + a) 
— — 


-l3/2 


(to) 


The quantity C is now evaluated. By equations (29) and (2l) 


= J(<P±,-8) dti'j = [y~ njg(9j) 


where 


(to) 


r P+“ -A [a ( x+y) 2+cy2j 

g(q>) = / J(<P,3) dti = K e L - 1 


ay 


(to) 
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by equation (39) . Thus, 

I Aac(Ap) 2 

8(<p) e ~ c+a {k5) 
As before, the sum (43) is split into two parts: 


C " Y~ 85 51 ^ n kS(9i)g(q>k) +51 n 2 f 2 (cpi) 

i^k i ' 


= 5l(t) sin *Pi sin ‘Pk s(<Pi)g(<Pk)A 2 + 

ij^k 



s AN sin cp^ 
2 


4 * 


.s AN sin 


(46) 


and taken "to the limit A— ^0 so that the second term of the second sum 
in equation (46) vanishes: 


C 



sin cp dtp 


(47) 


If Ns » 1, 
equation (47) 
by cos $ Q , 


C 


which is always implicitly assumed, the second term in 
becomes negligibly small. Replacing again sin cp 



1_ 

ac 


s cos t3q Kjt\ 2 
2A ) 


(48) 
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With this value and from equation (42): 


2 

5 = 


-u 

Ns cos -a \fsiAa. 2^^(a + cY’^c'^ 


(49) 


or, substituting the definitions ( 38 ), 


S 2 - 


( 4q sin^, -a + p ) 2 

2^‘TTs cos^ -S sin -8 \Ja$Tt ^4a sin^ -Q + 


(50) 


DISCUSSION OF RESULTS 


An immediate consequence of equation (50) is that the mean deviation 
" 2 ) 1/2 

i J is inversely proportional to the square root of the number N of 
diffracting crystal grains. This result was to he expected on general 
grounds . 


For a given number N, the mean deviation increases like 




or 


(equation 5) if a. is small. For small crystal dimensions p the 
deviation is inversely proportional to p. For large crystals, the mean 
deviation tends toward: 




1 


sin -5 



cos -d 


(51) 


/~2 \l/2 

Follow the change of ) as a function of a. If $ is large, 
sin $ x 1 may be assumed and equation (50) becomes: 


| 2 = 


Nsko2^ 2 x cos^ 3 sin $ 


+ P 


a 


(52) 


■OlH 
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The deviation decreases monotonically with increasing a. It will have 
twice its minimum value 



yl/2 


■a sin dy 


when 


or 



P = 


/ 3 

j8.l6 x 56 


For copper Kbc-i radiation, [3 = 0. V58 X 10 


-8 


so that 


P 


= 3 X 10 


-6 


centimeter 


(53) 


For values much smaller, the deviation increases rapidly; hut for larger 
values of p, the deviation decreases only slowly. While it is desirable 
to make the crystal size large for a given number H, an increase of the 

grain dimension beyond 10 ^ centimeter = l/lO p does not increase the 
accuracy noticeably . 

So far, the number of crystal grains has been considered as inde- 
pendent of the grain size. In most cases, however, the irradiated volume 
rather than the number N is constant. When the sample is sufficiently 
thick, this volume V is the product of irradiated area 0 and the 
effective depth of penetration t given by the absorption of X-rays. 

Then, ^ 

y P 3 W = V (5*0 


^Thls assumes that the grains occupy a continuous volume which is 
not literally possible for spherical grains. However, the assumption of 
spherical form is arbitrary and in bulk material the grains are actually 
shaped so as to fill a continuous volume . 
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or 


N 



(55) 


If sin -9 £ 1 is again put in the factors including a in equation (50), 


i 2 - 


a\Aa + 


cos 2 d sin d \/p" • sVjt 2 kQ 


(56) 


For a given irradiated volume, the mean deviation increases monotonically 
•with the grain size . In particular, for small grain sizes (4a « 3) , 

the mean deviation (s 2 )"*"^ is proportional to the grain radius. 3 

The quantity of main interest is the error in lattice-parameter 
measurement Ad/d. From Bragg's law 


Xq = 2d sin d Q 

the following relation is obtained by logarithmic differentiation 


a 


Ad 

d 


COt dg Ad 


(57) 


If -3 q is measured in the fashion described, that is, as the centroid 

of the intensity curve, the error A3 equals the quantity The mean 

deviation becomes: 



cot -9 



X 0 <x 4a sin 2 d + |3 

2 TA„ (6lna) 3/2( i|OBin 4 4 +p )3/4 


(58) 


^This result should not be taken literally for extremely small grain 


sizes (p sy 10 ^ cm) 
lines and fails when 


because our calculation supposes nonoverlapping 
the lines become too broad. 
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by equations (50) and (55)* An example of this relation is given in 
table I. Clearly, the relative error Ad/d again increases with 
increasing grain size and for small grains the error is proportional to 
the grain radius. 


How the relative error depends on the Bragg angle will now be 
considered. For very small angles d , the relative error becomes very 

la* 0 e because of the factor (sin -3 ) 3/^. However, for medium large angles 
(~8o°) there is only a slight variation of the error, and there is no 
noticeable advantage in using extremely large angles such as 87-1°. High 

angles merely serve to bring out more clearly the details of the erratic 
intensity fluctuation. 


and 


c 

If the grain radius is of the order of 10~° centimeter or larger 
•g is of the order of 45° or larger, then 4a sin^ 3 » 3 and 


Ad 

d 


' x o U/a a 3A 

£3/\(sin 3)5/2 


( 59 ) 


As an example, take the irradiated area to be 1 square millimeter, the 

effective depth of penetration as 10~3 centimeter, p = 5 x 10 - ^, and 
sin d ~ 1. Equation ( 59 ) reduces to 


Ad ^ 1 3/2 

T ~ 2 p 


or, for crystal radii 




This is an accuracy 


about 10 times better than that reported in experimental literature 
(reference 3). It maybe concluded that, by eliminating all sources of 
error not due to the grain statistics, the accuracy could be improved 

by a factor of 10 , for grains of the order of 10" ^ centimeter and 1- square- 
millimeter irradiated area. 


Only speculations can be made as to the main cause of error in 
actual measurement. It would seem that the most important factor not 
inherent in the crystal grain statistics is the graininess of the 
recording film. 


Armour Research Foundation 

Chicago, 111., Jan. 30, 1948 
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TABLE I 

VARIATION IN MEAN ERROR OF LATTICE-PARAMETER 
MEASUREMENT WITH GRAIN SIZE 

[iron-cobalt Kcx radiation; 6 , 80.64^; 0.8l x 10~"^ 

irradiated area, 1 sq mm; effective depth, 0.02 mm 
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Figure 1.- Schematic diagram showing diffraction by a single crystal grain. 



Figure 2.- Schematic diagram showing spatial orientation of a crystal grain 
by the position of a crystallographic direction on the unit sphere. 





